Atomic quantum gases in Kagome lattices 
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We demonstrate the possibility of creating and controlling an ideal and trimerized optical Kagome 
lattice, and study the low temperature physics of various atomic gases in such lattices. In the 
trimerized Kagome lattice, a Bose gas exhibits a Mott transition with fractional filling factors, 
whereas a spinless interacting Fermi gas at 2/3 filling behaves as a quantum magnet on a triangular 
lattice. Finally, a Fermi-Fermi mixture at half filling for both components represents a frustrated 
quantum antiferromagnet with a resonating-valence-bond ground state and quantum spin liquid 
behavior dominated by continuous spectrum of singlet and triplet excitations. We discuss the 
method of preparing and observing such quantum spin liquid employing molecular Bose condensates. 
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During the last 30 years, condensed matter physics has 
devoted a considerable interest to the issue of frustrated 
quantum antiferromagnets (QAF) (cf. 0,0)- In 1973 
P.W. Anderson proposed |3j the resonating-valence-bond 
(RVB) for the ground state of QAFs, where all spins are 
paired into singlets. RVB states exhibit neither the stan- 
dard antiferromagnetic Neel order, nor the spin-Peierls 
order (for which singlet pairs are ordered in space). Re- 
cent extensive numerical studies have shown that the 
RVB physics characterizes the spin 1/2 Heisenberg an- 
tiferromagnet on the 2D Kagome lattice in 2D (see Fig. 
la) 0,0. The spectrum of this system has a very pecu- 
liar structure: the energy gap between the ground state 
and the lowest triplet state, if any, is predicted to be 
very small (of order J/20, where J is the spin exchange 
coupling). This gap is filled with low- lying singlet states. 
Their number scales as 1.15^ with the number of lattice 
sites N. For temperatures above the triplet gap, the spin 
correlations decay rapidly in space, but have a very slow 
temporal behavior (s(x,0)s(x,t)) cx l/t os 

A very illuminating analytic insight into the physics 
of QAF in Kagome lattice has been recently obtained 
by Mila and Mambrini 0, who considered a trimerized 
Kagome lattice (TKL, see Fig. lb). Such a lattice is 
a triangular lattice of trimers with intra- (inter-)trimer 
couplings J and J' <C J, respectively. In Refs. a non- 
trivial mean field theory has been formulated that pre- 
dicts correctly the number, the form, and the spectrum 
of singlet excitations, which correspond to a restricted 
set of short-range RVB states. For J' <C J the theory 
predicts a triplet gap (2/3) J'. All these theoretical find- 
ings do not yet have a clear experimental confirmation 
in condensed matter systems, and it is thus desirable to 
seek for other possible testing grounds. 

Such novel testing grounds could be provided by the 
atomic physics of ultracold quantum lattice gases, which 
is one of the most rapidly developing fields of AMO 
physics nowadays. Following the proposal of Jaksch et 
al. 0, Greiner et al. were the first to demonstrate 
the superfluid-Mott insulator (MI) transition in a lattice 



Bose gas, predicted earlier in Ref. [l(|. Atomic physics 
and quantum optics offer in this context a new and very 
precise way of preparing, manipulating and detecting the 
system under investigation. 

In this Letter we show that using superlattices tech- 
niques 01 j it is possible to create a 2D optical trimerized 
Kagome lattice, and to control in real time the degree of 
trimerization (i.e. the ratio of J'/ J). The physics of 
cold atomic gases in such an optical lattice is described 
quite generally by various versions of the Hubbard model, 
and the energies and couplings defining such models can 
be calculated using solid-state methods (ti ght binding 
approximation, Wannier function expansion |12j). Such 
unprecedented possibility motivates us to consider three 
kinds of quantum gases in the TKL: i) A Bose gas with 
repulsive interactions; for this case we predict the ap- 
pearance of Mott-insulator phases with fractional fill- 
ings v = 1/3,2/3; ii) A spinless Fermi gas with nearest- 
neighbor interactions; such gas appears for instance natu- 
rally as a st rong -interaction limit of the Bose- Fermi Hub- 
bard model 13] . In this case fermions (or, more generally 
composite fermions consisting of bare fermions coupled to 
bosons or bosonic holes) attain boson mediated interac- 
tions. Interestingly, at 2/3 filling the Fermi gas behaves 
as a frustrated quantum magnet on a triangular lattice 
with couplings dependent on the direction of the bonds; 
iii) Finally, we consider a Fermi-Fermi mixture at full to- 
tal filling Ni = N2 = N/2, where Ni is the number of 
fermions of each species. Such a system in the strong 
coupling limit is equivalent to the spin 1/2 Heisenberg 
antiferromagnet in zero magnetic field, and the physics 
of this system is well-described by the RVB model. We 
discuss in details how to prepare a system in low-lying 
singlet states, and how to detect its properties. 

In the following, we consider the atoms in a 2D optical 
lattice in the xy plane, being strongly confined (magnet- 
ically or optically) in the z direction. In order to form a 
Kagome lattice one can use blue detuned lasers, so that 
the potential minima coincide with the laser intensity 
maxima. A perfect triangular lattice can be easily ere- 
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ated by two standing waves on the xy plane, cos 2 {ki^f) : 
with k\ t 2 = fc{l/2,±V3/2}, and an additional standing 
wave cos 2 (&3r + cj>), with = fc{0, 1}. The resulting tri- 
angles have a side of length 2^ j\[Zk. By varying <j> the 
third standing wave is shifted along the y axis, and, in 
principle, a Kagome pattern could be realized. 

However, this procedure presents two problems. First, 
3 lasers on a plane cannot have a mutually orthogonal 
polarization, and consequently undesired interferences 
between different standing waves occur. This can be 
avoided by randomizing the relative orientation of the 
polarization between different standing waves, or by in- 
troducing small frequency mismatches, which, however, 
have to be larger than any other relevant frequencies. 
Second, due to the diffraction limit, the ratio £ be- 
tween the separation between maxima and the half- width 
at half maximum (HWHM) is 4. Because of that, for 
any <j>, the three potential minima forming the Kagome 
triangles cannot be resolved. This can be avoided by 
using superlattices. For instance, one may substitute 
each standing wave (i — 1,2,3) by a laser potential 
(cos(fc.;r) + 2 cos(fciF/3)) 2 . For this potential £ = 7.6, and 
the generation of a perfect Kagome lattice for <j> = ir/2, 
and a modestly TKL for 57r/12 < <f> < ir/2 is possible. 
Using another superlattice: (cos(kif) + 2cos(kif/3) + 
4cos(fc.;r/9)) 2 , one can reach £ ~ 14, and generate for 
(f> = 7r/4 a strongly TKL (Fig. 1(b)). The necessary laser 
configuration is sketched in Fig. 1(c). 
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Figure 1: (a) Ideal Kagome lattice for <j> = n/2; (b) TKL using 
4> = 7r/4. This lattice can be generated using 3 standing waves 
with 7r/3 angle between themselves. Each standing wave is 
generated by 3 lasers with a configuration shown in (c); (d) 
enumeration of spins in a trimer and of neighboring trimers, 
used in Eq. (3). 

The proposed scheme allows one to realize the trimer- 



ized and the ideal Kagome lattice, and transform one into 
the other and back by changing 4>. One is thus tempted 
to study the physics of ultracold gases in such lattices. 
In the following, we consider three different physical 
systems: a Bose gas, a single component Fermi gas 
with nearest-neighbor interactions, and a two-component 
Fermi gas. The physics of these systems is quite gener- 
ally described by a Hubbard model [j|. The parameters 
of the model include intra- (inter-)trimer tunneling (hop- 
ping) J (J'), on-site interaction energies V for different 
species, and nearest- neighbor interactions U (£/'). These 
parameters can be controlled by modifying the laser de- 
tunings and intensities, as well as the configuration of the 
superlattices (the phase </>). 

To calculate the parameters J, J',V, . . . we determine 
first the Bloch functions for the problem. Since there are 
3 lattice sites per Bravais lattice, the lowest band splits 
actually into three different bands • For each of the 
bands we calculate the corresponding Wannier functions 
(localized at each one of the Kagome triangles) , and com- 
bine them then by means of an orthogonal transformation 
to construct functions localized at each lattice site. Once 
these functions are known, the tunneling and interaction 
coefficients are calculated as in Ref. 0- 

We consider first the Bose gas, which is described by 
the corresponding Bose-Hubbard Hamiltonian 
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where the tunneling and nearest-neighbor interaction co- 
efficients and Uij take here the values J and U for 
intra-, and J' and V for inter-trimer hopping, (i is the 
chemical potential, m = b\bi, and b\ are the annihila- 
tion and creation operators for bosons at the site i. In the 
limit when the boson number Nb < N, i.e for filling fac- 
tors v < 1, and for strong on-site interactions 7> J, J', 
only one particle may occupy a lattice site. 

In the strongly trimerized case (J',U' <C U < J) 
the system will enter a trimerized Mott phase with the 
ground state corresponding to the product over (inde- 
pendent) trimers. Depending on the particular value of 
p, = (/i — U)/(2J + U) we may have (p < -1), 1 
(-1 < p < 0), 2 (0 < p < 1) or 3 (1 < p) bosons 
per trimer, i.e. filling factors v = 0, 1/3, 2/3 or 1 bo- 
son per site. For fractional filling, the atoms within a 
trimer minimize the energy forming a, so-called, W-state 
[13: \W) = (jOOl) + |010) + |100»/V3 for v = 1/3, and 
\W) = (|110) + |101) + |011))/V§ forv = 2/3. It is worth 
noticing that VF-states themselves have interesting appli- 
cations for quantum information theory(c.f. jlfij). 

Generalizing the Landau mean-field theory of Ref. ^3 j 
we obtain the phase diagram in the J' = J' /(2J+U) and 
p plane with characteristic lobes describing the bound- 
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aries of the Mott phases, given by J' = — l)/2 for 
|A*I > 1, and J' = (3/2)[^|(l - \fl\)/(4 - \fl\) for | M | < 1. 
Observations of this Mott transition require tempera- 
tures T of the order of J', i.e. smaller than J and U. 
Assuming that U is of the order of few recoil energies , 
that implies T in the range of tens of nK. The results for 
J < U are qualitatively similar. 

The spinless Fermi gas in the trimerized Kagome lat- 
tice is described by the Fermi-Hubbard Hamiltonian 

H F h = - y^(Jfj///j+h.c.)+y^ UijUiUj-^fmi, (2) 

(ij) (ij) » 

where as before faj and Uij take the values J and U for 
intra-, and J' and V for inter-trimer hopping, /i is the 
chemical potential, n, = fa fa, and fa are the fermionic 
annihilation and creation operators. In the following, we 
enumerate the sites in each trimer 1, 2, 3 clockwise start- 
ing from the upper left site. We denote the 3 different 
intra-trimer modes by / = (fa + fa + fa)/V3 (zero mo- 
mentum mode), and f± — (fa + z±fa + z±fa)/V3 (left 
and right chirality modes), where z± — exp(±27r«/3). 
The intra-trimer Hamiltonian acquires the form Hi ntra — 
—3J fa f + U[(n — fa) 2 — fP]/2, where n is the total fermion 
number in the trimer, whereas fb = (fi — J +U /2)/U. In 
the strongly trimerized limit, we have (depending on the 
value of p): (/2 < 1/2 - 3J/U), 1 (1/2 - 3J/U < j± < 
3/2), 2 (3/2 < fj, < 5/2) or 3 (fi > 5/2) fermions per 
trimer. Obviously, the cases with and 3 fermions per 
trimer are not interesting. In the case of 1 fermion per 
trimer or less (filling factor < v < 1/3), at low temper- 
atures (T < J) the fermions will preferably occupy the 
lowest energy zero momentum mode / forming a gas of 
/-fermions on a triangular lattice. This gas will have a 
tunneling rate J'/3, and a coupling constant for nearest- 
neighbor interactions U'/9. Depending on the sign of 
U' we expect here similar behavior to that discussed in 
Ref. 0|, i.e. the appearance of a superfluid phase or 
fermionic domains for U' < 0, and Fermi liquid, or den- 
sity wave phases for U' > 0. 

The situation is more complex in the case of 2 fermions 
per trimer, since the second fermion may then have left, 
or right chirality, while the first fermion occupies the / 
state with certainty. The system becomes equivalent to 
a quantum magnet with interactions that depend on the 
bond directions, described by the Hamiltonian: 

2U' 6 

H^t = -5- XI s *(&->i) s J ( 3 ) 
i i=i 

where the nearest neighbors are enumerated clockwise 
as shown in Fig [Jd). In Eq. J2J) we have s((f>) = 
cos(ct))s x + sin(0)s. y , where s x = (/{/_ + /l/+)/2 and 
s y = - /l/+)/2. The angles cj> are <^_i = 

4>i^e = 0, 0.^2 = 0i^3 = 27r/3, </>i^4 = <j>i^5 = -2n/3, 
<j)i~,i = 4>2^i = -2-7T/3, (j)3^i = <j)A~,i = 0, and 4> 5 ^ t = 
(j)^ t = 27r/3. 



The couplings between the left and right modes are 
effectively ferromagnetic for U' > 0, and antiferromag- 
netic otherwise. Inter-trimer hopping J' contributes first 
in order (J') 2 /U, and is neglected here. For U' > the 
spins align ferromagnetically in the xy-planc in the clas- 
sical ground state. Quantum fluctuations around this 
state are represented in the semiclassical approximation 
by gapless magnons (cf. 0). For U' < the classical ap- 
proximation yields a Neel-type ground state with planar 
order in which on every triangle there is exactly one spin 
directed along the <f> = 0, 27r/3, — 27r/3 directions. Due to 
the anisotropy of the Hamiltonian there is no continuous 
rotational invariance in spin space. Therefore this planar 
ground state is unique (up to discrete 27r/3-rotations), 
and hence the semiclassical excitation spectrum in this 
case is gapped with a gap of order \U'\. 

The observation of the quantum phases in this model 
requires U' < J,U, and T < U'. If the fermionic inter- 
actions are due to dipolar forces 0|, then U may be of 
order of few recoils, and T in the range of 10-100 nK. If 
fermionic interactions are due to hopping induced effects 
in a Bose-Fermi mixture (as in Ref. 13]), then T in the 
range of 10 nK will be necessary. 

Finally, we consider a Fermi-Fermi mixture with half 
filling for each species, i.e. a spin 1/2 Hubbard model 

H FF = - J2 Jij (flfj + Pi fj + h-c.) + X Vn ^ ( 4 ) 

(ij) * 

where the operators fi and fj (fi and //) are the cre- 
ation and annihilation operators for the two components, 
Tii = fafi (n% = fifi), and, as above, J id = J (Jq) for 
intra- (inter-) trimer hopping. In the strong coupling 
limit, J , Jo <C V (t — J model) Q, Hff reduces to the 
Heisenberg antiferromagnet 

H = J ^2 Si ■ s.j + J' ^2 Si-sj, (5) 

<»,J>intro (*,j><r.ter 

where J = 4Jg/V, and J' = 4J'q/V, and s — (s x ,s y , s z ), 
with n — h — 2s z , fa f — s x + is y , and fa f = s x — is y . 

In the strongly trimerized limit 0, the total spin in 
the trimer takes the minimal value, i.e. 1/2, and there 
are four degenerate states having s z = ±1/2 and left or 
right chirality. The spectrum of the system in the singlet 
sector consists of a narrow band of low energy states of 
the width of order J', separated from the higher singlet 
(triplet) bands by a gap of order 3J/4 (2J'/3). 

Assuming that we can prepare the system in a sin- 
glet state at J' < T < J, then the density of states of 
the low lying singlet levels can be obtained by repeated 
measurements of the system energy. The latter can be 
achieved by simply releasing the lattice, so that after tak- 
ing care of the zero point energy, all of the interaction en- 
ergy is transformed into kinetic energy. In a similar way 
we can measure the mean value and the distribution of 



4 



any nearest-neighbor two-spin correlation functions. To 
this aim one has to apply at the moment of the trap re- 
lease a chosen nearest neighbor two-spin Hamiltonian and 
keep it acting during the cloud expansion (for details see 
|l7|). In a similar manner we can measure the spectrum 
of triplet excitation, by exciting a triplet state, which 
can be done by flipping one spin using a combination of 
superlattice methods and laser excitation 0] . The mea- 
surement of the singlet-triplet gap requires a resolution 
better than J'. 

A similar type of measurements can be performed in 
the ideal Kagome lattice, when J = J'. In this case, the 
singlet-triplet gap is filled with singlet excitations . By 
varying (j), one can transform adiabatically from strongly 
trimerized to ideal Kagome, for which the final value of 
J will be smaller than the initial J, but larger than the 
initial J'. In that case, the system should remain within 
the lowest set of 1.15^ states that originally formed the 
lowest singlet band. The singlet-triplet gap, if any, is 
estimated to be < J/20, and should be measurable using 
the methods described above. 

A possible way to prepare a singlet state in the TKL 
with T < 3J/4 could employ the recently obtained 
Bose-Einstein condensates of molecules consisting of two 
fermionic atoms a ^ temperatures of the order of 10 
nK. Such BECs should be loaded onto an ideal and weak 
Kagome lattice. Note that the molecules formed after 
sweeping across a Feshbach resonance, are in a singlet 
state of the pseudo-spin s. This can easily be seen, be- 
cause the two fermions enter the resonance from the s- 
wave scattering channel (i.e. in the symmetric state with 
respect to the spatial coordinates), and thus are in a 
singlet state of the pseudo-spin (i.e. antisymmetric state 
with respect to exchange of electronic and nuclear spins) . 
Since the interaction leading to the spin flipping at the 
Feshbach resonance |20j is symmetric under the simul- 
taneous interchange of both nuclear and electronic spin, 
then the formed molecule remains in a pseudo-spin sin- 
glet. The typical size of the molecule is of the order of 
the s-wave scattering leng th a, and thus can be modi- 
fied at the resonance |2l|, being chosen comparable to 
the lattice wavelength. Growing the lattice breaks the 
molecule into two separate fermionic atoms in neighbor- 
ing sites in the singlet pseudo-spin state. In this way, a 
singlet covering of the Kagome lattice may be achieved, 
allowing for the direct generation of a RVB state Q. 

Summarizing, we have shown that by employing cur- 
rently available superlattice techniques it is possible to 
create in a controlled way a trimerized Kagome lattice. 
An ultracold Bose gas in such lattice exhibits novel Mott 
insulator phases with fractional filling v = 1/3,2/3. A 
single-component Fermi gas with nearest neighbor inter- 
actions for v = 1/3 behaves as a Fermi gas in the under- 
lying triangular lattice, whereas for v = 2/3 it becomes 
a non-standard ferro- or antiferromagnet in the trian- 



gular lattice with direction dependent bonds. Finally, 
for a Fermi-Fermi mixture in a Kagome lattice, which 
is described by an antiferromagnetic Heisenberg model, 
we have shown the possibility to measure the spectral 
properties of the latter system. This opens the way of 
analyzing in a completely novel and fascinating setting 
one of the paradigmatic problems of condensed-matter 
physics, the physics of random valence bond frustrated 
antiferromagncts . 
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